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ABSTRACT 


We show that the Belinfante-Swihart (BS) theory can be 


reformulated in a representation in which uncharged matter 
responds to gravity in the same way as in metric theories. 

The BS gravitationally modified Maxwell equations can also 
be put into metric form to first order in the deviations of 
the physical metric from flat space, but not to second order; 
consequently the theory is nonmetric except in first order. 
We also show that the theory violates the high precision 
Eotvos-Dicke experiment, but cannot be ruled out by the 


gravitational precession of gyroscopes. 
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I. INTRODUCTION AND SUMMARY 


This paper analyzes the most complete and extensively developed non- 
metric theory that exists: the 1957 theory of F. J. Belinfante and 

12 3 

J. C. Swihart. ’ ’ Belinfante and Swihart constructed their theory as a 

k 

Lorentz symmetric linear field theory which would be easily quantized. 
However, as we shall show, in terms of measurable quantities the theory 

has all the nonlinearities of typical "curved-spacetime" theories. More- 

\ 

4 

over it is nearly a metric theory: We construct a new mathematical repre- 

sentation which has metric form to first order in deviations of the physical 
metric from flatness, but does not have metric form to higher orders. 

Section II gives a brief summary of the original BS representation. 
Included are discussions of nonlinearities and the behavior of rods and 
clocks. Section III presents our new mathematical representation of the 
theory. Section IV gives a prescription for obtaining the post- Newtonian 
Limit ’ of the theory and Sec. V considers various experimental tests. 
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Contrary to previous calculations it is found that both the geodetic and 

Q 

the Lens-Thirring precessions of gyroscopes cannot distinguish BS from 

General Relativity (for a particular choice of adjustable parameters). 

' 9 

However, using results of another paper, we show that the failure of the 
theory to be metric at second order causes a violation of the Eotvos-Dicke 
(ED)^ experimental results. Our calculations confirm the Belinfante- 
Swihart conclusion that their theory agrees with the three classical tests 
of gravitation theories, (perihelion shift of Mercury, bending of light by 
the sun, and redshift of light), and also agrees with the weak equivalence 
principle 11 (WEP) to first order. 
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II. THE BELINFANTE-SWIHART REPRESENTATION OF THEIR THEORY 


a. Lagrangian and Equations of Motion 


The original representation of the BS theory is Lagrangian based 11 but 

4 

is not in generally covariant form. In this section we generalize, in a 
trivial manner, the original representation so that it is generally covari- 
ant. The dynamical equations are obtained by extremization of the following 
action: 

1 = J + J ' (1) 
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Equations (1-5) describe the interactions of a collection of charged particles 
(labeled by A) with the electromagnetic and gravitational fields. Conven- 
tions and definitions for the above are the following: 

(i) We use units such that c = G = 1. 

(ii) TJog is a Riemann flat background metric (absolute gravitational 
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field ). In some coordinate system it therefore takes on 
Minkowski values, = diag. (-1, 1, 1, 1 ). All tensorial indices 
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occurring in Eqs. (1-5) are raised and lowered with q _. 

Gp 

(iii) Greek and Latin indices run through 0-3, 1-3, respectively, 

(iv) a, f, K are adjustable parameters. 

(v) h = h is a symmetric second rank dynamical gravitational 

)+ 

field. 

(vi) The world line of particle A is parametrized by an arbitrary, 

monotonic parameter A^ which varies from - « to + oo. Particle 

A is described by its coordinate z^ and its "velocity and 

momentum variables" a^ 1 and which are all functions of A^. 

(vii) The electromagnetic field is described by the tensor fields 

A and H = - H . 
p pv vp 

(viii) T^ v is a "stress-energy tensor" for particles and electro- 
magnetic . fields. (The bar above is used to distinguish it 
from a different "stress-energy tensor" definded in Sec. IV.) 
(ix) Slashes denote covariant derivatives with respect to the flat 


background metric t) . 
(x) q ■= determinant of q . 


Equations (5a) and (5b) are decomposition equations for T^ and b^. The 

dynamical variables which one varies independently in the action are h (x), 

p v 

a A^(\)> n A^^A^ A ,,( x )> and H 1M1 (*)- Variation of the matter 


pv' 


variables yields the following dynamical laws 
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[BS, I, (210] 
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tBS, II, (10)) 
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b. Nonlinearities in the Theory 


Linear gravitational field equations do not preclude a nonlinear form 

. I 

for the response of particles to gravity. The BS theory is an example: 

Eqs. (16) and (17) endow the canonical variables and with gravita- 
tional contributions. Consequently the equation of motion for a particle, 

Eq. (10), is nonlinear in the gravitational field h . Indeed, although 

p v 

the BS theory is often called a "linear" theory, its "linear" first-order 
matter Lagrangian produces qualitatively many of the nonlinear effects of 
General Relativity (GRT) for example (see Secs. Ill and IV). Hence one 
should be cautious in the labelling of theories as "linear" or "nonlinear" 
on the mere basis of the linear forms of their gravitational equations. 


c. Behavior of Rods and Clocks 
3 

In the third paper of their series Belinfante and Swihart quantize 
the theory and obtain a gravitationally modified Dirac theory. We remind 
the reader that all nonmetric theories must exhibit explicitly the manner 
in which all the laws of physics are changed in the presence of gravity. 



Belinfante and Swihart find that, in the case of a static spherically 
symmetric (SSS) gravitational source, the standard solutions to the unmodi 
fied Dirac equation are related to those in the presence of gravity in the 
following way: ^ 

cp o (x o , t) = Ncp(x, t) 


x = C x 
~o 


[ BS, III, (78)] 
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Here the subscripted quantities are those in the absence of gravity, cp is 
the electron wave function, U is the Newtonian gravitational potential for 
an SSS source and "a" is the previously mentioned adjustable parameter. 

The coordinate system is one in which = diag. (- 1, 1, 1, l). The 
energy eigenvalues, i.e., E in cp(x, t) = cp(x) exp(- iEt/h), are shifted in 
the presence of gravity: 

E Q = (1 + U)E [BS,III, (82)] 

— a result following essentially from Eq. (l4). It is Eq. (16) which pro- 
duces qualitatively fhe correct redshift. Equations (12) and (13) also 
indicate the effect of gravity on the coordinate sizes of atoms. Consider 
the expectation value of the coordinate size of an atom: 

< r > = J* l«p(x» fc ) | 2r <* 3 x • 

Using Eqs. (12) and (13) we obtain 
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According to Eq. (16), the coordinate ticking rate of an atomic clock 
decreases in a gravitational field: u> = (l - u)a) Q . According to Eq. (IQ) 
the coordinate size of a rod made of atoms increases in a gravitational 
field: l = [l + U(| a' Q * Since a » to agree with the light bend- 

ing experiment (see later sections) the above results are the same, to 
first order in U, as one obtains in GRT, using an "isotropic, post- Newtonian" 
coordinate system. ^ 


III. ATTEMPTS TO PUT THE THEORY INTO METRIC FORM 

The BS theory is a Lagrangian-based, relativistic theory of gravity.^ 
Therefore, according to a theorem proved in Ref. 4, it is a metric theory 
if and only if the "nongravitational part" of its Lagrangian; 

^NG = ^M + ^1 

can be put into universally coupled form. ^ Let us try to achieve universal 
coupling by a change of variables, i.e., by introducing a new mathematical 
representation of the theory. 

a. Particle Part of Lagrangian 

Begin with the terms in ^ NG that refer only to particles and define 
the following tensors: 

V * », M - s K 

= A/ ^ = T 


a 


( 19 ) 

( 20 ) 


Then, from Eq. (7), obtain the relation 


a v = Z^ V dz^/dA 


( 21 ) 
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Equation (21), which is obtained after variation of the Lagrangian, suggests 
that one define a new variable v* 1 to replace a^ 1 in the Lagrangian: 


v _ A v , Lt 
a = A v 
. H 


( 22 ) 


Then, the relation v^ = dz^/d A will presumably turn out to be an Euler- 
Lagrange equation. Using Eqs. (l9)-(2l), bring the particle portion of 
the Lagrangian into the form 
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In obtaining Eq. (24) from Eq. (23) we have performed the integrations over 
d^x and thus all of the spacetime functions should be evaluated at the 
particle position z^. 

If we now define an "effective metric" 

Sop ' < l - Kh > 2 W - V 1 - 2KH) + "oS + ° (h2) 


(25) 


Eq. (24) takes the universally coupled form, with g^ being the only gravi- 
tational field occurring in I .. . Variation of n then yields the desired 

part. p 


relation 


= dz^/dA 


(26) 


To make our results look simpler, we expliticly introduce Eq. (26) 
into Eq. (24), thus eliminating II completely and obtaining 
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Variation of Eq. (27) yields equations of motion which, by the use of 
Eqs. (6), (19)- (21), can be shown to be identical to the BS equations of 
motion, Eqs. (10). Equation (27) is the familiar "metric theory" action 
principle describing the interaction of charged particles with the gravi- 
tational field g^ and the electromagnetic field A . 

b. Electromagnetic Part of Lagrangian 

It will now be shown that, to first order in h , the electromagnetic 
Lagrangian can also be put into metric form. Change variables from H 

pv 

to an antisymmetric tensor F by 

V” v (1 + l h + ¥ h2 ) <» < h > 

- 2F a( - 2 f a/[A] A + 0(f1,3 > • < 28 > 

Equation (28) is simply the result of an inversion of Eq. (8). Square 

brackets around indicies denote antisymmetrization of indices (with the 

usual normalization of a factor of -i). Variation of F in. the new 

2' pv 

Lagrangian presumably will yield the relation 


F = A . - A . 

pv v |p p I V 


(29) 


Substitution of Eq. (28) into the electromagnetic portion of the action yields 
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where 


i- 1 * ^ \ n^d - \ h> + 

If one now uses the inverse of Eq. (25), i.e., 

g ce , ^08 _ h os + 2Kh 0S + 0(h2) 

and 

(- g) 1 / 2 = (- n) 1 / 2 [1 + h(| - UK)] + 0(h 2 ) 
one finds Eq. (32) can be written as 
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Thus has universally coupled form at 0(h); at 0(h ) deviations occur, 
arising from the L corr term in Eq. (33). Variation of F^ v in Eq. (33) 
yields the desired relation between F^ v and A^, i.e., Eq. (29). Completely 
equivalent equations are obtained if Eq. (29) is now substituted into 
Eq. (33), yielding 


- - f 18 ”)' 1 . 

r vvv g0<5gllV( - 8)1/2 dk * + w. 

(35a) 

= - ( 16it)~ 1 ( 

r v^t- 8)1/2 dU * + l cor R . • 

(35b) 


The relation given in Eq. (29) is now understood to hold in Eqs. (35). 
Since we now have constructed a second metric g^ (the "physical metric"), 
indices on all quantities except the constituents of g^ (q^, h op> £^p) 
henceforth will be raised and lowered with g^. Equation (35), aside from 
the 0(h ) correction term, is recognized as the electromagnetic Lagrangian 
for metric theories. Thus the BS theory is a metric theory at first order 
but nonmetric at all higher orders (in h). 
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c. Summary of Our New Representation 

Our new representation of the BS theory is summarized succinctly in 

Table I. In particular one sees that for uncharged particles the theory 

is metric to all orders in h, with playing the role of the "physical" 

k 

metric. When electromagnetic phenomena are included, and when one goes 

beyond first order in h, the theory is nonmetric (cf. in Table i). 

LURK 


IV. THE POST- NEW TONI AN LIMIT OF THE THEORY 

We now proceed to calculate the post-Newtonian (PN) limit of the theory. 
The PN limit is a perturbation solution of the gravitational field equations 
— expanding in the small quantities .occurring in the solar system, e.g. , 

= (macroscopic velocities .of bodies)^ = O(e^) 

11 1 

i 2 

U = Newtonian gravitational potential = 0(e ) 

2 

p/p = (pressure)/(proper density of rest mass) = 0(e ) 

2 

II = (internal energy density)/(rest mass density) = 0(e ) 

We refer the reader to Ref. 5 for further details of the expansion scheme. 

a. The Metric-Theory Approximation 
Using Table I, we write the field equations as 
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where we have used the usual definition (as in metric theories) 

ixv _ 2 ^Metric 

7 072 6g 

(- g) ' HV 


( 37 ) 


To PN order, the first term on the right-hand side of Eq. (36) is of order 

2 

First Term = (Total energy density) • e 
while the second term is of order (see Table i) 

2 

Second Term ce (Electromagnetic energy density) • e 

Since the electromagnetic energy of a substance is typically smaller than 

—3 

the total mass-energy by a factor < 10 , the second source term in Eq. (36) 

can be neglected at PN order, by comparison with the first. Similarly, one 

I I I ' ! 

can make a metric- theory approximation for the response of matter to gravity. 


For metric-theory (i.e., universally-coupled ) Lagrangians one always has 

1^ V = 0 

;v 

when the matter field equations are satisfied, where the semicolan denotes 

covariant differentiation with respect to the physical metric g . In the 

Qp 


(38) 


BS case 


T** V « h h 


Sh 


(39) 


pv 


so again one can conclude that effects resulting from the deviation in the mat- 

-3 

ter response equation from Eq. (38) will be 5 10 of PN effects. Thus for all 
PN phenomena we can neglect and treat the BS theory as a metric theory. 


b. From Point Particles to Perfect Fluid 


In one of their original papers Belinfante and Swihart, when solving 
their gravitational field equations with the sun as the external source, 
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use an ad hoc perfect-fluid stress energy- tensor for ? iV , rather than the 
expression given in Eq. (5). Their T is precise enough to yield an ade- 
quate treatment of the "three classical gravitation tests" but not precise 
enough to adequately handle such effects as the effective gravitational 
mass of gravitational energy (cf., "Nordtvedt effect" in Ref. 5). To 
avoid such problems, and to ensure self-consistency of the theory when 
dealing with gravitating sources in the solar system, we will build up the 
fluid BS stress energy tensor, T^ V , as an average over charged point par- 
ticles and their electromagnetic fields [cf., Eq. (27) and Table I]. 

The kinetic- theory procedure for constructing a perfect fluid out of 

interacting particles is the same in any metric theory as in general 

relativity — and the same in general relativity as in special relativity 

13 

("equivalence principle"). By following that standard procedure and by 
neglecting the resulting non-perfect fluid terms, we obtain the standard 
stress energy tensor: 

Tl = (e + p) tfu + pg^ 

Here u* 1 is a suitable macroscopic average of the microscopic particel 
velocities; e is the density of total mass-energy (rest mass plus kinetic 
energy of particles plus electromagnetic energy) as measured in the macro- 
scopic rest frame; and p is the similarly measured averaged pressure. 

c. The Parametrized Post-Newtonian (PPN) Formalism 

References 5 and 6 present a "parametrized post-Newtonian formalism" 
in which the PN limit of every metric theory is summarized by the coeffi- 
cients of various integral functions in its metric. These coefficients, 
the so called PPN parameters, are obtained by the previously mentioned 



perturbation solution (PN limit) of the gravitational field equations. We 
have constructed such a solution for our new mathematical representation 
of the BS theory, using Eqs. (36), (40), and Table I. The details are 
spelled out in Ref. l4. (Actually Ref. l4 is the presentation of an exact 
gravitation theory closely related to the metric- theory approximation of 
the BS theory. ) We refer the reader to Ref. 14 and here quote only the 
PPN parameters of the BS theory. 

7 = 7 + 0(w) £ 2 = ° = 0(w) 

p = p + o(w) £ 3 = ° a 2 = 0(w) 

q = 0 = o 0^ = 0 . ( 4 i) 

Here 7 and P are given implicitly in terms of a and f by 1 

a = 1/(2 7+2) (42) 

f = 10 F + 6 7 P ■ 7 7^ - 87-6 
2(7 + 1 ) (3 7 + 5 - 4 p ) 2 


and to obtain the correct Newtonian limit, one must require 

16K 2 a - 4aK + a + 5f _ 0 

a(a + 4f) ~ * (44) 

By 0(w), we denote terms involving the cosmological boundary values of h^^ 

(see Ref. 14 for further details). Imposing Eq. (44) reduces the number of 
arbitrary parameters to two (a and f for example); so we may regard 7 and 
P as being arbitrary. For comparison, General Relativity has no arbitrary 
parameters and its only nonzero parameters are 7 = P = 1. 


i 
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V. EXPERIMENTAL CONSEQUENCES AND TESTS OF THE THEORY 


g 

In his 1972 Varenna Lectures, Will summarizes, within the PPN frame- 
work, the constraints which may be placed on a metric theory's parameters 
by current solar system gravitation experiments. As has been indicated in 

the previous section, the difference between the BS theory and a metric 

3 

theory for PPN- type experiments is less than one part in 10 . For most 
experiments the microscopic internal energies play a minor role — e.g., 
it is the macroscopic rotation of the earth which produces the macroscopic 
Lens-Thirring precession of gyroscopes. For such experiments the BS theory 
is effectively a metric theory to a much higher accuracy than indicated 

above. In summary, so far as PN experiments are concerned, to the precision 

! I i 

of the technology of the 1970' s, the BS theory is accurately summarized by 
the values of its PPN parameters, Eqs. (4l). We refer the reader to Ref. 8 
for the experimental consequences of those values. Here we merely point 
out a few salient features. 

Perhaps the most important feature is this: If the 0(w) terms in the 

parameters are sufficiently small, and if the arbitrary parameters are chosen 
so that 7 = F = 1, then the PN predictions of the metric- theory approximation 
to BS are the same as the PN predictions of general relativity. In parti- 
cular, the predictions for the "three classical tests" are the same as 
Belinfante and Swihart* themselves deduced by complicated calculations. 

a. Preferred Frame-Effects 

For the coordinate system in which q is Minkowskian, it is natural to 
set the boundary values of h to zero when treating the solar system, as was 
done originally by Belinfante and Swihart. However, the correct way to 
determine the boundary values of h is through the solution of the cosmological 
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problem. If the solution produces nonzero cosmological boundary values of 
h, then those values will effect certain of the PPN parameters [cf., 0(w) 
terms in Eqs. (4l)]. In the case of the BS theory the presence of such 

I 

terms is a direct consequence of the presence of the "absolute gravitational 

4 8 

field" q (cf., Table i), and leads to various preferred- frame effects 

such as anamolous solid earth tides and contributions to the perihelion 
shift of mercury. We refer the reader to Ref. l4 for a more complete dis- 
cussion of the derivation of such effects in the BS theory. 

b. Precession of Gyroscopes 


We specifically mention this experimental test only because there seems 
15 

to be some confusion as to the prediction of the BS theory. Using for- 

I i : ' N 

mulas from Ref. 8 and the BS PPN parameters, Eq. (41), one ohtains for the 

precession of the spin S of a gyroscope orbiting the earth 


dS 

. dir s*§ 

where 

~ ” -LENS - THIRRING ^ -GEODETIC 
0 — 7 + 4 + 0(w)][.05" of arc/year] 

L- To t 

1 v 

n =— {1 + 27 + 0(w)J[7" of arc/year] 

G 5 


( 45 ) 


(46a) 

(46b) 

(46c) 


Thus the results of the upcoming (to be launched before 1977) Everitt- 
16 

Fairbank gyroscope experiment can only place upper limits on the cos- 
mological boundary volues of h [cf., 0(w) terms in Eqs. (46)] for a 

v 

given choice of 7 . 
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c. The Weak Equivalence Principle and Eotvos-Dicke Type Experiments 

We conclude by considering the Eotvos-Dicke (ED) type experiments, ^ 

which test gravity so precisely that they fall outside of the PN realm of 

precision. Braginsky's^ recent version Of the ED experiment shows that 

the difference in accelerations of test bodies of aluminum and platinum in 

12 

the gravitational field of the sun is smaller than one part in 10 . Such 

1+ 

a result represents a strong validation of the weak equivalence principle 

2 2 

(WEP). Consider the contribution of electromagnetic energy at order F h 
(see bottom of Table i) to the gravitational mass and acceleration a of a 
test body: 


- n 


Ji I /electromagnetic energy ! .2 
^ total energy 


E.M. energy 
^ total mass 


(47) 


2 2 = 

where h « U and g — VU. For platinum, the following relation holds: 


(E.M. energy/total mass) 10 


and the Newtonian potential due to the sun at the earth is 


U « 10 


■8 i 


Equation (47) and the above numerical estimates indicate that the ED experi- 
ment can distinguish between the BS theory and its metric- theory approxima- 
tion (cf., * n Table i). All metric theories satisfy WEP identically. 

The BS theory, however, as is shown in Ref. 9, predicts 


^Pt " ^Al 

»» 

g 

in clear violation of the Dicke 
The reader is referred to Ref. 


6 X 10 
10 


-11 


( E.M. energy \ y 
total mass J 


m 


ii 


and Braginsky versions of the experiment. 
9 for complete details as to the derivation of 
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Eq. (U8) from considerations of particles interacting with gravity and 
elec tromagne t ism. 
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TABLE I 


A New Mathematical Representation of the Belinfante-Swihart Theory 


1. Gravitational fields: 

a. Absolute field q 

b. Dynamical symmetric second rank tensor h 

c. "Physical" metric g 


2. Nongravitational variables: 

a. Particle coordinates z. 

~A 

b. Electromagnetic vector potential A 

c. Affine parameter of particle world lines 

I 

3. Gravitational field equations: j 

a. Flatness of q: Riemann (q) = 0 

b. Field equations for h obtained by variation of in Lagrangian 

below , 

2 u 

c. Decomposition equation for g: = (l - Kh) where we have 

defined (5^ T - ^ , K is an arbitrary constant, h = q^h^, 

and indices are raised and lowered on h^, with q^. 

4. Influence of gravity on matter: 

Equations for A, obtained by variation of those quantities in 

Lagrangian 

5. Lagrangian density: 
a * ^ = ^G + ^NG 

■>• *0 - - O 6 ")' 1 ( ah '‘ a |aVa |a + £h |a h ' a ) ^ 


18 



TABLE I (continued) 


• ^NG ’ ? S 


dz“ dz /' 1 / 2 
A A 


dz. 


“ m A\ g Qp dA A dA A 


+ e A A 


A V dA 


A J 


d \5 (? ~ 5 A ) 


( 16lt ) V g ) + *^corr. Metric + ^corr. 


where , the "correction term" in the Lagrangian, which repre- 
sents the amount by which the purely electromagnetic portion of the 
Lagrangian fails to have metric form, satisfies 

=Wr. = o^ 2 ) t 8ee E( i s * W] 
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